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The electron dephasing time in a diffusive quantum dot is calculated by considering the interac- 
tion between the electron and dynamical defects, modelled as two-level system. Using the standard 
tunneling model of glasses, we obtain a linear temperature dependence of 1/t^, consistent with the 
experimental observation. However, we find that, in order to obtain dephasing times on the order 
of nanoseconds, the number of two-level defects needs to be substantially larger than the typical 
concentration in glasses. We also find a finite system-size dependence of t^,, which can be used to 
probe the effectiveness of surface-aggregated defects. 

PACS numbers: 73.23.-b, 72.70.+m, 73.20.Fz 

I. INTRODUCTION 



Interference of the electron's paths in a mesoscopic 
conductor results in various quantum phenomena such 
as the universal conductance fluctuation, persistent cur- 
rent, and weak localization. In all these phenomena, 
the dephasing time appears as a typical time scale 
over which the electronic trajectories have interference; 
weak localization correction to conductivity — for exam- 
ple, is conventionally used fotj the experimental determi- 
nation of the dephasing timeH. In the moderate temper- 
ature range, experimentally-determineco values of in 
diffusive metals are found to be in excellent agreement 
with the theoretical predictions of due to electron- 
electron interactionEl. While it is theoretically expected 
that Tcf, — > oo as T — > in the absence of other ex- 
ternal sources of dephasing, is found tOpSaturate at 
low temperatures in almost all experimentst _ including 
the recent carefully-performed experimentsolJ. This se- 
vere discrepancy between theory and experimental ob- 
servation of low temperatuna|-saturation has fast be- 
come a topic of contiupversyffl surrounding the ques- 
tion whether the ideaQ and the theoryu of zero-point 
fluctuations of the electromagnetic field created by the 
electron-electron interaction as a source of dephasing 
are tenable on general grounds. This poses a serious 
problem as zero temperature dephasing of electrons has 
been argued to be relevant to the problems of per- 
sistent current in normal metalslij, the low tempera- 
ture metal/insulator, quantum-Hall/insulator and su- 
perconductor/insulator transitionaljLJ, and transport 
through r garipus normal- metal/superconductor hybrid 
junctionsoE3; but the most unsettling consequence is 
the negation of the fundamental premise upon which the 
theories — and hence our understanding — of metals and 
insulators are based: the many-body Fermi liquid pic- 
ture. 

Among various sustained efforts to find a zero temper- 
ature dephasing mechanism other than electron-electron 
interaction, dephasing due to dynamical defects inside 



the conductor has been recently arguecOliU to be impor- 
tant to the saturation problem. Low-energy excitations 
of the dynamical defects are usually modelled by two- 
level systems (TLS). Invoked some three, decades ago, 
first byp4nderson, Halpcrin and VarmaEa, and also by 
Phillipsli3, the tunneling model of TLS has been quite 
successful in explaining various anomalies in the acous- 
tic, dielectric and thermodynamic properties of struc- 
tural glasses and other amorphous solidsElx 

Imry, Fukuyama and Schwabli3 have recently suggested 
that the saturation behavior may have the same origin as 
the 1// conductance noise, arising from the two-level de- 
fects. Zawadowski, von Delft and RalphEJ have argued 
that the apparent saturation of may be caused by 
the two-channel Kondo effect due to reiectron-TLS scat- 
terings. However, it was pointed outturn that hysteresis 
or switching behavior, expected from the effects of TLS, 
was not observed in experiments. In addition, various 
concentration-dependent Kondo-like bulk trends antici- 
pated in these theories were also not observed in the ex- 
periments. 

In this paper, we investigate the role of two- level de- 
fects in the dephasing of electrons in quantum dots. In 
the recent experiments from the Marcus group, Huibers 
and coworkers have observed_the saturation of dephasing 
time in open quantum dotsEil below 0.1 K along with a 
strong temperature dependence above 0.1 K. In addition 
to this experiment, saturation of in qmmtum dots has 
also been reported in other experimentsEfcj, although 
the physical meaning of in these set-ups are not the 
same. If two-level defects are responsible iqr the satura- 
tion of dephasing time in this experimentsEil, it is natural 
to suppose that just above 0.1 K the linear temperature 
dependence should be explained by two-level defects as 
well. Our calculations indeed show that the dephasing 
rate due to the TLS does have a linear temperature de- 
pendence. However, we find that the magnitude of the 
dephasing by two-level defects is too small to explain the 
experimentally observed dephasing time of nanoseconds. 
This implies that either other mechanisms are more effec- 



1 



tive or surface-aggregated two-level defects play a dom- 
inant role; defects on a disordered surface are likely to 
have unusual distributions in their splitting energies; we 
suggest that the surface defects can be experimentally 
probed by measuring the size dependence of dephasing 
time. 





FIG. 1. Two types of electron interferometry involving 
quantum dots (a) Both paths are inside the dot. (b) One 
of the interferring path is outside the dot. 

Consideration of two-level defects in quantum dots for 
their dominant role in dephasing is motivated by the 
experimentally observed tell-tale signs of TLS in quan- 
tum dots: hysteresis and switching behaviors, which 
have been iu |Jact detected in various quantum-dot 
experimentsEScH unlike in the experiments on higher di- 
mensional diffusive metals. In a quantum dot, usually 
the Thouless energy Et is theJargest energy scale un- 
like the diffusive metallic caselm. Therefore, the results 
obtained in the diffusive metallic casesEm cannot be ap- 
plied (even after the appropriate dimensional considera- 
tions) to quantum dots. Thouless energy is defined by 
Et = hD I L 2 , where D is the diffusion constant and L is 
the typical system size. 

Dephasing generally describes the loss of coherence or 
suppression of interference. Accordingly, it is important 
to know which kind of paths are considered before defin- 
ing the typical time scale of the loss of interference along 
these paths. In Fig. 1, we show two kinds of electron in- 
terferometers involving quantum dots. In this paper, we 
are concerned with pairs of time-reversed paths which 
return to the origin in a diffusive system (Fig. 1. (a)). 
These time-reversed paths enclose magnetic flux; their 
interference manifests in the weak localization correc- 
tion to conductivity. These paths are chosen, for the 
problem at hand, because their contribution to conduc- 
tivity does not vanish even after disorder averaging. In 



the interferometry discussed in this paper, change of the 
mean conductance at a finite field from its zero field 
value, Sg = ((g)) B ^ — ((5))s=o' can be used to ex- 
tract the dephasing time from experimental data, where 
((...)) means disorder averaging. Using the phcnomelog- 
ical random matrix theory (RMT) in Refs.oE^I, can 
be defined, for instance, by the formula 



Sg: 




(1) 



where N is the number of channels connected to a quan- 
tum dot. Although, the formulae for the conductance 
change are model dependentBEj, the difference in the 
equations in these models are not significant to the inter- 
pretation of measured in experiments^. In this work, 
we will—refer to which were obtained through Sg as 
in RefEj without discussing how Sg is related to any 
further. 

Our calculation of dephasing time is similar to the 
calculation of dephasing time by Stern, Aharonov and 
ImryEj. Based on the interference of two time-reversed 
trajectories, we calculate a typical time scale over which 
the environmental state remains in the initial state. De- 
phasing rate due to electron-electron interactions in 
diffusive quantum dots has been calculated by Sivan and 
coworkersElt 
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Si 



Et) ' 



e » k B T; 



(2) 



where e is the excitation energy of the particle and Si is 
the mean level spacing. However, it should be noted thai 
the direct application of Eq.(|J) to experimental dataEJ 
is difficult, because it is not meaningful to estimate the 
temperature dependence of t$ by merely replacing e with 
~k B T in Eq.0B. 

While the interference paths considered here are 
inside the rdfft-i as depicted in Fig.l.(a), there are 
experiments£-3H associated with other types of interfer- 
ence involving one path through the dot and another out- 
side the dot (Fig. 1. (b)). Interference appears as a peri- 
odic Aharonov-Bohm oscillation rather than the fluctua- 
tions in the electric current as in the case of Fig.l. (a). 
In experiments, dephasing manifests in the nonpideality ( 
< 1) of the visibility of the interference patterns. Theo- 
retically, dephasing rate involved in the geometry 
in Fig. 1. (b) is calculated-as the level broadening 
In these considerationscili23, level broadening is due to 
the gauge fluctuations of a OD-quantum state where the 
electron dynamics inside the dot is neglected. There are 
other ways to d&fine dephasing time. In the work of Sivan 
and coworkers^!, dephasing time is considered to be the 
quasiparticle life time, measured in the tunneling experi- 
ments. In the work of Bird and coworkersca, conductance 
fluctuation was used to measure at very strong mag- 
netic fields. Nevertheless, we will refrain from discussing 
these different views of dephasing time in this paper. 
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The organization of the paper is as follows. In section 
II, we describe how dephasing time is calculated in a gen- 
eral framework. In section III, interaction between the 
two-level defects and electrons is discussed. In section 
IV, the return probability is calculated in order to ob- 
tain Trf,. In section V, we discuss dephasing by two-level 
defects with widely distributed energies. In section VI, 
dephasing by two-level defects with narrowly distributed 
energies are discussed. 



II. INTERFERENCE AND DEPHASING OF 
PARTICLE'S TRAJECTORY 

Let us consider the event that the electron is at the 
position r at an initial time t = 0, and it arrives at the 
position r' by diffusive motion after a time To- The en- 
vironmental state changes from r\ to rj' in this process; 
the corresponding probability amplitude of the event is 
P(r',r, r/, 77; t ). 

The description of the suppression of interference in 
electron's paths by the electron-TLS interaction can be 
considered in two different approaches: 

(i) The electron in the two different paths produces 
two different time-dependent electric fields on TLS, 
thereby TLS go to different states, which sup- 
presses interference. 

(ii) The fluctuating dipole moment of TLS produces 
the time-dependent electric field, thereby the elec- 
tron in the two different paths gains random phase, 
which also suppresses interference. 

In general, these two approaches are not equivalent, 
because the presence of the electron induces a back re- 
action from the TLS environment. However, in the pres- 
ence of weak interaction between the particle and the 
environment, it is known that either the two description 
are equivalent, or at least they give the same dmhas- 
ing rate up to the second order in the interactioncl In 
this work, we use approach (i). Following the scheme of 
Chakravarty and Schmidts, we use semiclassical approxi- 
mation on particle's trajectory and we consider quantum 
mechanical evolution of the TLS (environment) states. 
We further assume that the TLS environment does not 
influence the classical paths of the electron, therefore the 
diffusive electron motion comes from only static disorder. 
Under certain conditions, the two-level defects might be 
able to effectively change the semiclassical paths of the 
electrons. In that case, one may estimate by calculat- 
ing electron-TLS inelastic scattering time. However, 
begins to lose its meaning as a dephasing time, since we 
lose the semiclassical picture of the electron's path. 

The tunneling motion in the TLS environment is to be 
described in a fully quantum mechanical way. To this 
end, we consider the time-dependent potential V(r(t)) 
exerted by the moving electron of the path r(t) on a two- 
level defect. The probability amplitude is given by 



Kr'.r.^^ro) = ^A j (r' ,r-,T )e iS * ( V '\Uj(T )\r,), (3) 
3 

where Aj and Sj are the corresponding amplitude and 
action of a classical electron's trajectory labeled by j. 
Uj{jo) is a time-evolution operator (in the interaction 
picture) of the environmental state associated with the 
electron trajectory Yj(t) 



Uj(ro) = Texp 



V!(Tj{t),t)to 



where T is the time-ordering operator and 



V I {r j {t),t)=e* tl ™ t V{r j {t))e-* 



rH env t 



(4) 



(5) 



The probability P(r', r, 77; To) of finding the particle at 
r' after time To, initially at r with the environment in 
the initial state |m, is given by the sum of the absolute 
square of the probability amplitudes over the final states 
of the environment; 



P v (r',r;T ) = J dr]'\p{r' ,r,rf ,rj;T )f 
= ^|^(r',r;T )| 2 



(6) 



+ J2^Ale^- s ^(r 1 \ut(To)U J (To)\r 1 ) 1 (7) 

using the completeness relation for the environmental 
states. 

The return probability P rj?) (To) of an electron is defined 
by the probability of finding the electron at position r 
after time to which was initially at the same position 
with the environmental state 1 77) ; 



PrM=P ri (r',r;T )\ T ,^ i 
= J2\A J (r,v 1 r )\ 2 



(8) 
(9) 



I- ir- LJ -\ r i) + 



where j T denotes the time- reversed path of j. The first 
term in Eq. ([)]) is termed as the classical return probabil- 
ity -P r dass (To). The second term comes from the interfer- 
ence of the pair of time-reversed paths. The remaining 
terms that do not appear in the above equation vanish 
upon ensemble averaging over disorder due to the random 
differences in its classical action Sj — Sk for k ^ j,j T - 
The coherent part P^Jtq) of the return probability is 
the second term in Eq. (j9J) ; 

p rt fa) = E l^(r,r ; To)| 2 Re<7 ? |C/;(To)C/ J T(To)|7 7 ). (10) 



Now, the dephasing time can be defined as a time scale 
for P coH (tq) to vanish with a decreasing function such as 
exp (— to/t^). But for the present purpose, the particular 
exponential form of time dependence is not needed. 
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III. ELECTRO N-TLS INTERACTIONS 

We consider two-level tunneling systems (TLS)0 as 
the environment for an electron in the quantum dot. Let 
us first consider TLSs which have asymmetry energy A 
and the tunnel splitting energy A . The TLSs are as- 
sumed to be randomly distributed over the dot with their 
electric dipole moments randomly oriented. We will as- 
sume the dipole moment is not too strong so that we 
do not have to consider interaction among the TLS. The 
density of the TLS will be assumed to be not too high 
so that multiple scattering events between the electron 
and the dipoles can be neglected. Within this approxi- 
mation, we calculate the return probability of an electron 
in the presence of a single TLS, thereby we extend the 
reults to the case of many randomly distributed TLSs. 
The Hamiltonian of the TLS can be written in terms of 
the localized wave functions of the double well potential 
and also in terms of the eigenenergy basis 



Htls = 





' A 






' E 






v A 


-A 




K 


-e) 



(11) 



where E = A 2 + Aq and the transformation denoted 
by the arrow means localized wavefunction representa- 
tion — > eigen wavefunction representation. The dipole 
strength operator p is defined in the eigen wavefunction 
representation: 



P = Po 



1 

-1 



Po 



A/E A /E\ 
A /E -A/E J 



(12) 



where po is the dipole moment when the particle is lo- 
cated in one of the wells of the defect potential. In 
the following sections, we will use the eigen wavefunc- 
tion representation in which Htls is diagonal. The TLS 
Hamiltonian will be used for the environment Hamilto- 
nian H env = H TLS . 

TLS dipole at the position R feels the electric field 
E(R) produced by the moving electron. The result- 
ing interaction energy can be expressed by the operator 

V(r 3 (t)) = -pn ■ E(R) = -pn ■ V R y c (|R - T d (t)\), (13) 

where p is the dipole moment operator for the TLS, which 
is along the direction of unit vector n. V c is Coulomb 
interaction potential; 



V c (K-r 3 (t)) = 



e*|R-r,-(i)| IA^ q 



(14) 



where e is the electric charge, e* is the dielectric constant 
of the dot material, L is the linear system size, and d is 
the spatial dimensionality of the dot (d=2, 3). Here, v q 
is the finite Fourier transform of Coulomb potential; 



dre- lqr — i-r 
dot e *l r l 

2ne , , , 

4-7re , , , 



(15) 
(16) 
(17) 



We use the discrete values of q = ^-(m,n 7 k) for the 
three-dimensional case, and q = ^(m, n) for the two- 
dimensional case where Lm, and n are integers. By in- 
serting Eq.Q into Eq.(|||), we get 

V(rj(t)) = -^P^2 ( in ' q « q e~ 4q ' R )e iq ^ (t) . (18) 
q 

Derivation of the Eq.(|i~8|) is based on the semiclassi- 
cal approximation of the electron's motion and the un- 
screened dipole moment of TLS. From a pure quantum- 
mechanical point of view, one can consider the TLS- 
electrOjtL. interaction similar to its treatment in metalic 
glasses^. When the TLS has two positions R ± = 
R± d/2, the purely quantum mechanical TLS-electron 
interactions Vqm are written ascil 



Vqm = Td - ]T 2i M ,e- iq R sin(q • d/2) 

X Y, C l C k+q, 



(19) 



where ct (c^) is the electron creation(annihilation) op- 
erator with momentum k, and /i q is the Fourier trans- 
form of the ionic potential. However, the pure quantum- 
mechanical approach is not reliable, because the concept 
of dephasing becomes ambiguous as we leave the semi- 
classical approximation, which has been pointed out in 
Ref.H. Here, we merely get some useful informations 
by comparing the quantum mechanical Hamiltonian in 
Eq.(|l9|) to our semiclassical potential in Eq. (|l8|) . Since 
q • d < < 1 in quantum dots (this is true-generally speak- 
ing, when the Fermi wavelength is much larger than d), 
sin(q • d/2) » q • d/2 = n • q|a. Eq.(fL9|) can be under- 
stood as the interaction between the TLS-dipolc coupled 
to an effective electric field produced by the electron. 

One of the two differences between Eq. (|lj) and Eq. (|l9| ) 
is that the electron interacts with the ion through a 
screened interaction \i q in Eq.(|l^) rather than the di- 
rect Coulomb interaction v q as in Eq.(18). The specific 
form of /i q is not known though it is expected to be less 
than Vq. If a screened interaction [i q is used, then the 
calculated dephasing rate l/r^ would be smaller than 
that with the unscreened interaction v q . In this work, we 
will use Vg instead of \i q . The second difference, which 
is rather important, is the quantum mechanical nature 
of the electron motion in Eq.(19). V(rj(t)) in Eq.(|l§|) is 
understood as (ipj(t, t)\VQM\ipj(r,t)) , where tfjj(r,t) is 
the time-dependent wavefunction which describes a wave- 
packet corresponding to the trajectory j; 
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k 



(20) 



At a finite temperature T, the wave-packet state 
|^> 3 -(r,t)) will consist of mostly the eigenenergy states 
with energies limited within tp ± fc^T. Therefore, the 
time dependence of e iq r ^*- ) will be limited by the fre- 
quency windows \w\ < ksT. Though we use the semi- 
classical approach in Eq.(p"g|), the frequency cut-off will 
be performed in our semiclassical calculation later. 



IV. THE RETURN PROBABILITY 

For the simplicity of calculation, we suppose that the 
TLS is initially in the state |— >, an eigenstate of the 
Hamiltonian in Eq.(|lT|) with the eigenenergy EL (the 
case with |+ > can be calculated in a similar way). The 
time evolution operator Uj(t) (Uj(t)\ — ) = c + (t) \ + ) + 
c_ (t) | — )) of the TLS corresponding to the electron path 
j can be written as 



c+(f) 
c_(i) 



Texp 



dfVi(rj(t)) 



(21) 



where 



Vi(rj(t)) = 





Vj( rj (t'),t') 




Vj{ rj {t'),t') 




(1- 


Vi(r 3 {t'),t>) 









(22) 



To find the corresponding time evolution of the TLS 
states for the time-reversed paths, (Ujt (t) I — ) = d + (t) | + 
) + <i_ (t) | — \), one obtains a similar form by using 
r jT (t) = Vj(T -t). Expanding (- \u] {t )U jT (tq)\ -) = 
c* + {To)d+{ro) + c*L(to)cI~(to) up to the second order in 
interaction V, we get 



Re(-|C/t(r )f/^(r )|-) 



1 + 

1 r° , 



dt 1 [cos(fi(* + t' - t )) - cos(0(i - t'))] 
x(-|V(r J -(t))|+><+|V(r i (f ))|->. (23) 



where Ufl = E = E + - EL = ^ IS. 2 + A§, and we used 
the relation ^(rjfiVjt = T/(i\-(f)) in the last equality. 

Now from Eq.(18), Eq.(|23|) and Eq. (|Io|) , one can get 
the coherent return probability in 



ErVo) = £l^(r,r;T ) 



3h 2 L 2d 



Here, we omitted the subscript rj in P£ °" for the TLS 
state, because both of the initial states of the TLS \rj) = 
| ± ) give rise to the same expression. P^ oH (tq) in Eq.(|2J) 
is the value averaged over the TLS position. Here, we 
used the disorder average over TLS << exp[iR • (q + 
q')] >>= <5 q ,-q'- The factor 3 in Eq.(p4[) comes from the 
average over the orientation of the TLS dipoles. 

In the case of ^/Dtq > I (I is the mean free path), the 
sum over the classical paths which appears in Eq.(|2~i|) 
can be jiraittep as a path integral using the Wiener 
measurecBBEj. The path integral can be calculated as 



[cos(0(i + t' - t )) - cos(fi(t - t'))] 



,(*')) 



^|^(r,r;r )| 2 e 



iq-(rj(t)-r;,(i')) 



x(r )=r 



£>[x(r) 



(25) 



x(0)=r 



- £ exp(-£| q fT ) e -^N 2 - 2 ^' 



)\t-t'\ 



V 

|q'|<7r/; 

where we used the boundary conditions of a quantum dot 
by a rectangular box of the volume V = L x Lx (L or 
a << L for 2D dots). Several remarks are in order. Eqs. 
( p5| ) are valid only in diffusive regimes. The summation 
over the discretized momentum variables q and q' is un- 
derstood to be limited by it /I. The contributions from 
ballistic regime, which are supposed to be small when 
*/ 1 Dt^ >> I, are neglected in this work. The last equal- 
ity in the Eqs. ( p5| ) is obtained for r = 0. To simplify 
the calculation, from here on, let us consider the return 
probability of the particle at the origin r = 0. 

Then, the classical return probability P class (To) for r = 
is given by 

P class (r )=Y,\M^ro)\' 

3 

= 1 E exp(-£>|q| 2 T ). 



(26) 



|q|<f/2 



Inserting Eq. (^5|) into Eq. (|24|) , we get the coherent part 
of the return probability 



P^(r )^cxp(-J}|ql 2 r ) 
q' 







\( + \p\- 





dV 



t+ 



dt / dt' 
J a 



3H 2 L 2d 

[ cos fi(t+ - r ) - cos fir] e -^(lql 2 -2q-q')l* 



dt- (27) 



where we used the change of variables of t + = t + t' and 
t~ =t- t'. 

(24) We restrict to the ergodic regime r > t^, therefore 
significant contribution comes only from q' = 0. Then, 
we obtain 



5 



P co \t ) 



l r 
V 



1 



cos f2(t 



3?i 2 L 2d 

Tq) — COS ilt~ 



dt~ 



L + 



1 

2^ 



exp(iw|i |) 



ito + Dq 2 



(28) 



The frequency of the time-dependent electric field pro- 
duced by the electron is not infinitely large but has a 
cut-off. By assuming the electron to be in equilibrium 
with other electrons at temperature T, the high frequency 
cut-off of lo is given by k B T (|oj| < k B T); this is true 
at temperatures that are not too low. Note that be- 
cause of the finite size of the system, qv q = for q = 0. 
Therefore, there is no divergence at low frequencies and 
the low-frequency cut-off of lo does not play an impor- 
tant role. By integrating Eq.(^8|) with the condition of 
< fcflT, we obtain the coherent part of the return 
probability i 5 ™ 1 *, which decays as oc f — tq/t^ + ■ ■ ■. We 
have now defined as the dephasing time. The dephas- 
ing rate from a randomly distributed TLS with an 
asymmetry energy A and a tunnel splitting Ao is given 
by 



1 



2|< 



\P\ 



7*(A,A ) 



3U 2 L 2d 
2p 2 A 2 



E2 2 
q v, 



Dq 2 



q 9, 2 + {Dq 2 ) 2 ' ^ 



3h 2 D(A 2 + A 2 ) 



E 



2,1 



(30) 



0<|q|<7r/Z 



We used flru << f in the second equality in Eq.(30). 

obtained above is valid only when k B T > Ml 1 while 
in the other case P coh (To) is an oscillating function of r 
with a small amplitude. Therefore l/r^ (k B T < htt) is 
negligible. 

By inserting Eq.@ and Eq.@ into Eq.@ when 
k B T > ^/A 2 + A 2 , we get 



2p 2 e 2 ~ d A 2 



r (A, A ) 3h 2 De* 2 L 2 (A 2 + A 2 ) 



(31) 



where d = 2, 3 is the spatial dimension of the quantum 
dot, and 



s»= E 

0<m 2 +n 2 <(L/l) 2 
f 



1 



"3 — 



7T- 



E 



0<m 2 +n 2 +k 2 <(L/l) 2 



k 2 ' 



(32) 
(33) 



V. DEPHASING BY TWO-LEVEL DEFECTS 
WITH WIDELY DISTRIBUTED ENERGIES 

We can generalize 1 /t# to the case where the TLSs are 
distributed with a distribution function /(A, A ), 



- = V [dA [ dA 1 

Tcf, J J i>(A,A ) 



/(A, A c 



(34) 



By inserting Eq.(|31j) into the above equation, we get 
f _ L d - 2 2p 2 e 2 ~ d a z - d 



T 4> 



D 



where 



S(T) = J dA J dA c 



37Te 



A 2 

A 2 + A 2 



-S(T), 



(35) 



/(A, A a )6(k B T — y A 2 + A 2 ,), (36) 

where a is the thickness of the dot in case d = 2. It is 
interesting to note that at d = 2, the dephasing time does 
not depend on the dot area. 

To calculate r^, we use,|fche standard tunneling model 
for the two-level defectsliata. The essential postulate in 
this theory is the uniform distribution of the tunneling 
parameter A associated with the tunnel splitting Ao oc 
e~ A . The energy distribution function /(A,Ao) in this 
case is written as 



/(A,A ) = £. 



(37) 



Furthermore, it is also assumed that Ao has a nonzero 
minimum value Ao, m m- By applying this distribution, 
we find; 



S(T) = P 



k B T 



= 2FA , r 



f ^/(k B T) 2 -Al 

dA / dA 

J-^/(k B T) 2 -A 2 g 

F{k B T/ A ^ min ), 



A 



A 2 + A 2 



where 



Hz) 



dxt'ATl 



(38) 
(39) 



(40) 



Here, the above expression is valid for k B T < Ao^ max , 
which is an realistic and common assumption for the tem- 
perature below IK. Note that F(z) ~ z In z when z » 1, 
therefore in the case of k B T » Ao, m m, we expect the 
following temperature dependence 



1 ti <JzbT\ 

— Tln ( A h 



(41) 



which is closer to ~ T rather than ~ T 2 . 

Now let us estimate quantitatively. We. consider 
the experiments by Huibers and coworkersEIJ on two- 
dimensional ballistic semiconductor, quantum dots. The 
quantum dots in the experimentsEj are in the ballistic 
regime, while our is for diffusive quantum dots. How- 
ever, since the dephasing time is in the ergodic regime 
(t^ > to): the results for diffusive dots should be ap- 
plicable to the chaotic quantum dots in the ballistic 



G 



regime. The diffusion coefficient is obtained through the 
ergodic time scale and D ~ (ETh/^i)(^/2m*), where 
to* is the effective mass of the electron ( fop-GaAs, 
to* = 0.067to 6 ), and E Th /h ~ 30 for RefEI For 
ballistic dots, the Thouless energy is given by Tivp/L. 
For GaAs, h 2 e*/m*e 2 ~ IOtito. A reasonable size of the 
dipole moment is po ~ e x 10~ 10 to. The thickness of 
the two-dimensional quantum dot is roughly a ~ IOtito. 
By putting together S2 and ln(fcsT/Ao irn i„), which are 
roughly ~ 1 — 10, into Eq.(^), we find 

l/r - (10~ 16 - 10- 15 )m 3 s- 1 Pfe B T. (42) 

In order to obtain ~ Ins near T = 0.1K, the average 
concentration should be P ~ (10 48 - 10 49 )J" 1 to~ 3 . Al- 
though this number is not completely unreasonable, it is 
too large to be expected from, well-textured semiconduc- 
tors used in the experimentsEil. For comparision, we note 
that glassy materials possess a typical TLS concentration 
of P~ 10 45 - lO^J- 1 ™- 3 . 

One may anticipate a different temperature depen- 
dence which might show the saturation of by consider- 
ing the dissipative two-level system due to TLS-phonon 
interactions or incoherent two-level systems due to TLS- 
TLS interactions. However, it is very difficult to expect 
that the dephasing rate is enhanced by several orders of 
magnitude by such interactions. 

The large magnitude of P may be possible if a large 
enough number of two-level defects aggregate on the sur- 
face of the quantum dots. This possibility can be ex- 
perimentally checked by varying the system size and the 
dimensionality. Using our results, 

1/t* oc L d - 2 a 3 - d . (43) 

For example, for a 2D quantum dot, the dephasing rate 
I/t^ by "intrinsic" two- level defects will increase as the 
thickness a of the dot increases, whereas it will decrease 
with a for surface defects. 



VI. DEPHASING BY TWO-LEVEL DEFECTS 
WITH A NARROW ENERGY DISTRIBUTION 

Low-energy excitations exist in semiconductor crystals 
due to the tunneling of impurity ions between equivalent 
interstitial lattice sites. Due to the crystal fields, def- 
inite positions are preferred and a wide distribution of 
excitation energies is not expected; in glasses, the wide 
distribution arises because of structural disorder. How- 
ever, defects on the surface may result in a wider distri- 
bution of energies because of surface roughness. A single 
tunnel-splitting energy implies a narrow distribution of 
relaxation times such that the standard tunneling model, 
applicable to structural glasses, as discussed in the pre- 
vious section, is not validlfJ. 

In this section, we consider a well defined tunnel- 
splitting energy Ap rather than a wide distribution. The 



asymmetry may be uniformly distributed with a gaus- 
sian width Ai, usually determined from the experimental 
data. The distribution function is defined as 

/(A) = n TLS -^e- A2 / A K (44) 

n-TLS is the TLS density. 

The function S(T) defined in Eq.(|3^) in the expression 
for the dephasing rate 1/r^ is simplified to 

S(T) = J dA-g£2— f(A)0(k B T-y/A* + A*). (45) 

Note that the variable Ao is not integrated over, in con- 
trast to the case for the standard tunneling model; and 
the final result depends on Ao. Evaluation of the above 
integral yields 

S(T)~n TLS (k B T»A »A 1 ) (46) 
~ T~ n TLS {k B T » Ai » A ). (47) 

If temperature is larger than the energy scales of TLS, 
then it is possible to obtain saturation or temperature- 
independent dephasing rate l/r^. In realistic systems, 
Ao is usually a small fraction of Ai; thus the experimen- 
tally relevant limit is the second case, Ai 3> Ao, in ex- 
pression (|47|) . The dephasing rate can now be obtained: 

- ~ (10- 16 - lO-^riTLS^-mh- 1 . (48) 

T<f, Ai 

For to be on the order of 1 ns, the two-level defect 
density should be 

"TLS~^(10 24 - 10 25 )m- 3 . (49) 

Now let us estimate utls for a typical single-crystal 
system. Single crystal silicon structures have been stud- 
ied in this context in the temperature range of inter- 
est, below 1 K down to 5 mK. Both acoustic dissipation 
and heat capacity measurements on silicon resonators by 
Kleiman, Agnolet and, BishopCj (see the corresponding 
estimates by PhillipsEil and KeyescZl) find that the TLS 
density, utls ~ 10 23 to~ 3 , with an estimated value of 
Ai/Aq ~ 100. Now using the same value for Ai/A 
in the expression (^9|), the order-of-magnitude estimate 
of TLS density is found to be n TL s ~ 10 26 - 10 27 to" 3 . 
The required density needs to be at least three orders of 
magnitude higher than the typical concentration in the 
silicon structures to result in a TLS-induced dephasing 
time 70 ~ 1 ns. This is an unreasonably large number, 
even for the typicaLintentionally-doped semiconducting 
structures of silicone! Though, experimental studies of 
acoustic and thermal properties of gallium arsenide struc- 
tures/heterostructures for the effects of two- level systems 
have not been-done, in the temperature range of interest 
for dephasingeju, recent-studies on semi-insulating gal- 
lium arsenide resonatorsc3 suggest that the typical TLS 
density is comparable to that in silicon. 
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VII. CONCLUSION 

We have calculated the dephasing time by assuming 
the presence of two-level defects inside diffusive quan- 
tum dots. The temperature dependence of l/r^ is found 
to be roughly ~ T for widely distributed two-level de- 
fects in the standard tunneling model. We find that to 
explain the size of the experimentally-observed dephas- 
ing times, we need a large number of two-level defects. 
This number is substantially larger than that found in 
glassy materials (almost by three orders of magnitude). 
Therefore, it is hard to believe that the electron dephas- 
ing is dominated by the intrinsic two-level defects at low- 
temperatures. We have also calculated t$ from a distri- 
bution of narrow energy two-level defects, and we find 
a regime of temperature independent t^. However, the 
required number of two-level defects is too large as in the 
case of widely distributed TLS. The system size depen- 
dence obtained in our calculation can be used to check 
the possibility of surface defects which are probably ef- 
fective. Because of the large surface-to-volume ratio in 
quantum dots, it may be reasonable to assume that most 
of the defects are surface-aggregated. It will be interest- 
ing to estimate P or utls required for the observed low- 
temperature charge noises of quantum dots and compare 
to the values from dephasing time. Unfortunately, we are 
not aware of any quantitative theory for the quantum-dot 
charge noise arising from the two-level defects. 
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Note added.- In a recent paper, Aleiner, Altshulcr, and 
GalperinEJ have analyzed the relevance of TLS for elec- 
tron dephasing. Although they use a different approach 
and evaluate for different systems (metals not quan- 
tum dots), their conclusions are similar to ours — that is, a 
substantially large concentration of TLS, P, much larger 
than the typical values in metallic glasses is required for 
the quantitative explanation of,-the saturation observed 
in experiments on metallic wiresta by two-level systemslla. 
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